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ABSTRACT 


“a This report describes a Fortran IV program capable of obtaining a 
finite element numerical solution to the free axisymmetric vibrations 


of a cylindrical type shell of arbitrary shape. 
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INTRODUCTION 


Present flextensional designs may be grouped into one of five different 
classes (1). In order to derive a mathematical model for any of the above 
mentioned flextensional transducers, it is essential to determine the 

# vibrational characteristics of the shell itself, A finite difference approach 
for each shell's dynamic response has been investigated at the Center for 
Acoustical Sutides in the Department of Mechanical and Aerospace Engineeering 
at North Carolina State University. An attempt in improving the widely 
used finite difference method led to the present investigation of the 
finite element method which is widely used in de’%sermining stresses and strains 
in Aerospace Research and Development. 

The program utilizes a finite element analytical model to calculate the 
total Mass and Stiffness matrices of an arbitrary shape cylindrical type 
sheli (2). .After the total Mass and Stiffness matrices are calculated, 
displacement boundary conditions are introduced by removing appropriate 
rows and columns corresponding to points of the shell which are rigidly 
restrained from motion. Finally, subroutine NROCT is used for the numerical 
solution of the shell's free axisymmetric vibration matrix equation, and 


tne eigenvalues, frequency coefficients, natural frequencies, and eigenvectors 


are then printed out. 
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PROGRAM ULTLIZATION 


The program is written in Fortran IV (G or H level) to be used on an 
IBM 370/165 computer. The storage capacity required depends upon the 
number of finite elements desired. As an example, the sample program 
presented in this paper having a total number of 6 elements requires 


a storage of 100 K for the evaluation of the stiffness and mass matrices 


of each element. These numbers are stored on a disk. In addition, a 
a storage of 160 K, together with a time parameter of 4 minutes is 
i required for the main program for a complete run. A typical element 


& geometry and coordinate notation is given in figure 1. A geometrical 


— 
ete 


gi element assembly is given in figure 2. 
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Figure 2. Finite element assembly of a shell. 


Displacement vector 2: per node is given as follows: 


35 = My UL was us Das Qed Bu 9 
QH “Ox Ae” WY Is DY 


The program presented in Appendix A was written for the Class IV 
flextensional transducer shell, but it can be easily charged so it can 
be used for the case of any other similar type shell. Refore using 
this particular program, it must be noted that all element stiffness 
and mass matrices are evaluated outside the main program. Then tney 
ere put on a disk from whict ‘he main program picks them up for total 
mass and sviffness formulation and final dynamic solution. It is 


neceseury to use proper JCL cards. The program instruction in Appendix A 


MARE ie thee Rb ee SOE ROR ED EDS SE TH RISE Te ORIN a SR Rt ra ETE 


EVIE A 


must be read very carefully before running the program. A set of data cards 


is needed for the element notation and final mess and stiffness formulatico:. 


The program presented in Appendix A was written for figure 1. element 


notation and figure 2. element combination. All elements to form data 


for set 1. must be numbered according to figure 1 and figure 2. 
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m pe ~ * oe Sa 
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ELEMENT NO. 1 5 12 
ELEMENT NO. 4 9 5 6 10 
ELEMENT NO. 2 6 2 
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3 
ELEMENT NO. 5 106 7 
ELEMENT NO. 3 T3483 

8 


ELEMENT NO. 6 pe 4 12 


The boundary conditions are: Wy. =O b= 4) + 3,6,9, 12 
Veo > u 
Vy =O be " 
Ws =o be QlO ab ard 
Ww =O t= Wy 
Wy: = vey 


These are applied to the corresponding vector index notation 


WW eee Wy Vy Wy 
Node 1 1 ©) 3 * © 6 7 8 @ 
q y Node 2 10 11 212 #13 1 6015 16 «#17 «(18 
Node 3 19 20 21 22 23 2k 25 26 27 
Nede 4 2 € 30 31 @ 33 3 35 @ 

: Node 5 37 @ 39 40 @M ke 43 bb 

Node 6 46 47 48 49 50 SL 53 


Node 7 55 56 57 58 59 60 
Node 8 6h 67 69 
86 


® @ 


Node 9 ©) 75 76 
Node 10 33 8485 
Node 11 #2 93 9 95 @) 
Node 12 00) G63) 102 103 (0) 
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Therefore Set 3, correspdonding to the total number of the above mentioned : 


30 boundary conditions (NBC#30), is easily formed as follows: 


2-5-9 - 20 - 32 - 36 - 38 - WL 
45 - 65- 68- 72 - 73 - 74 - 77 - 78 set 9 i 
- 79- 81- 82- 87 - 68 - 91 - % - 97 ‘ 


100-101-104-105= 106- 108 


In order to form Set 2, the vector index notation is used again, but the 


Teta ne ON nae 


naw 


circled numbers that correspond to the boundary conditions are omitted. 


Therefore, the vector index notation is as follows 
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Node 1 BOYS ope ME) FO os, 
Node 2 @8 9 10 12 12 13 W 15 3 
Node 3 1617 18 19 20 21 22 23 2k ; 
Node 4 25 _ @ 27 _ @8) 29 30 : 
Node 5 GD)_ 62) 33 _ G35 3% _ : 
Node 6 67) 38 39 ho 4a 42 43 bh bs 
Node 7 4647 48 4 50 52 52 53 5h 
Node 8 55_ 6) 57 _ (8) 59 60 _ 

Node 9 eae cee Cs) 
Node 10 _@65 6 _ _ _ @6 

° Node 11 [5OU? BRB . . 
Node 12 oie 7) 09) 


8 

‘ o a 
3 Then Set 2 is formed as follows: 
: Ba hei? 606s B64 Biorgpe Bh 
a 37- 56- 58- 61- 61- 61- 63- 63 Q 
| ; SET i 
63- 6h- 64- 68- 68- 70- 7h- 7H 

a 76- 76- 78- 78- 78- 79 

3 Then the date cards consist of Set 1- Set 2- Set 1- Set 3 in that order. 

Furthermore if a coordinate transformation is desired, then matrix AKL (I,J), $ 

. 
: card MAIN 85, should be assigend properly and also cards DISK 27 and MAIN } 
ad 266, 267, 268 should be removed from the program. In addition, card MAIN 54 
a should be equal to 1. 
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FINAL REMARKS 


A 5 point numberical integration is carried out for each element's 
stiffness and mass matrix by using abscissas and weight functions of 


Legendre polynomials. After combining all the individual stiffness 


tae 


and méss matrices and introducing boundary conditions, the equation 


eS 


had 


of motion for the shell is written as follows: 
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Total mass matrix 
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st Nodal vector 
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Density 
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Modulus of Elasticity 
; FOISS= Poisson's Ration 


‘ Since subruutine NROOT is used to determine + >» the largest » is 


determined first, and as a xesult, the lowest natural frequency is 


calculated first 91a so is the corresponding eigenvector. 
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